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ABSTRACT 

We carry out a direct search for bar-like non-circular flows in intermediate-inclination, 
gas-rich disk galaxies with a range of morphological types and photometric bar classi¬ 
fications from the first data release (DRl) of the CALIFA survey. We use the DiskFit 
algorithm to apply rotation only and bisymmetric flow models to H a velocity fields 
for 49/100 CALIFA DRl systems that meet our selection criteria. We find satisfactory 
fits for a final sample of 37 systems. DiskFit is sensitive to the radial or tangential 
components of a bar-like flow with amplitudes greater than 15kms“^ across at least 
two independent radial bins in the fit, or ~ 2.25 kpc at the characteristic final sam¬ 
ple distance of ^ 75Mpc. The velocity fields of 25/37 (67.6t8;5%) galaxies are best 
characterized by pure rotation, although only 17/25 (68.0/(iob^) them have suf- 
hcient H a emission near the galaxy centre to afford a search for non-circular flows. 
We detect non-circular flows in the remaining 12/37 galaxies. We con¬ 

clude that the non-circular flows detected in 11/12 {91-Tlilfg%) systems stem from 
bars. Galaxies with intermediate (AB) bars are largely undetected, and our detection 
thresholds therefore represent upper limits to the amplitude of the non-circular flows 
therein. We find 2/23 (8.7/(2'g%) galaxies that show non-circular motions consistent 
with a bar-like flow, yet no photometric bar is evident. This suggests that in ^ 10% 
of galaxies either the existence of a bar may be missed completely in photometry or 
other processes may drive bar-like flows and thus secular galaxy evolution. 

Key words: 

galaxies: spiral, galaxies: kinematics and dynamics; galaxies: structure, surveys 


1 INTRODUCTION 

Current estimates indicate that 30% of nearby spiral 
galaxies are strongly barred in optical light, a num¬ 
ber which rises to ~ 50% if weak bars are included 


* E-mail: Lindsay.Holmes@rmc.ca 


llSellwood fc Wilk insonI 19931: Ba razza. J osee fc Marinoval 


200^ _ Aguerri. MendezAbreu fc Corsinil 20091 : 

Masters et al.l 1201 ll ). Bars are even more prominent in 


NIR imag es where the measu r ed ba r fraction is more than 
70% llMulchaey fc Regan _ 1997: lEskridge et al.l l2000l : 


IWhvte et al.ll20oT Sheth et al. 1201 ll ). The increase in bar 
fraction at NIR wavelengths has been attributed to a higher 
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prevalence of weak (AB) bars, which a re obscured by dust 
and star forminR regio ns in the optical llAthanassoulal[l992l : 
iMarinova et al.l[2009al l . Bars are therefore an essential struc¬ 
tural com£onent_whose_£ropertes drive secular evolution 
(e.g. lKormendv fc Kennicuttll20o3 . , and references therein). 
The presence of a bar within a galaxy allows for the redistri- 
bu tion of both the material and angular mom entum within 
it llAthanassoula, Machado fc Rodionov 2013li. Bars have a 


significant influence on gas flow ( Weinberg fc Kat3l2002ll . 

they are thought to b e an efficient mechanism for moving 
gas to central regions ([S hlos man. Frank fc BegelmanI 1 19891 : 
IShlosman, Begelman fc Fran^l990ll adding mass to th e 
bulge or creating secondary bars ( Friedli fc Martinetlll993ll . 
The incidence of barred galaxies is therefore important 
in understanding the secular processes that affect the 
basic structure of a galaxy and slowly change it over time 
dKormendv fc Bended 

Information on the morphological, photometric, 
and kinematic properties of bars has been obtained 
from many observational studies r anging from sin¬ 
gle objects to large surveys (e .g. Sheth et al.l 200^ 

jMarinova et al.l 


Aguerri, Mendez-Abreu fc Corsinil I200'. .... 

2009bl : ISimard et al.ll201ll : lAguerri et al.ll2015ll . Because of 


the relative ease of obtaining images for large samples of 
nearby galaxies, most observational studies of bars have 
focused on their photometric properties. While photometric 
studies of bars provide important insight into the role that 
bars play in galaxy evolution, the relationship between 
the properties derived from the light associated with bars 
and the dynamical impact that these bars have on their 
host systems is not well understood. In addition, bars in 
different dynamical states - young, flat bars and older bars 
that have buckled out of the disk, for example - exhibit 
only subtle differences in the ir photometric properties 
(e.g. &win_J^D_ebattist3 l2013l : iLaurikainen et akl l2014l : 
iMendez-Abreu et al.ll20l4) . 

A more direct measure of the dynamical importance 
of a bar in a disk galaxy is to measure the associated gas 
flows. However, studies of th e kinematics of bar-like flows 
have been performed (e.g. IWeiner. Sellwood fc William^ 
200ll: Hemandez_etalJ_ 20(35l : ISoekkens fc Sellwood 200?! : 


Sellwood fc Sanched 201Clh . and no thorough observational 


investigation of the correspondence between the photomet¬ 
ric properties of a bar and their kinematic signatures exists. 

The publicly released DiskFi10 software is specifically 
designed to fit for bisymmetric distortions in disk galaxy 
gas velocity fields. DiskFit fits non-parametric models to 
images or t o vel ocity fields as ori ginally described by 
iReese et al.l (120 071) for images, and Soekkens fc SellwoodI 
(l2007tl and ISellwood fc Sanched (l2010tl for velocity fields. 
DiskFit can fit both axisymmetric and non-axisymmetric 
models of disk galaxies, it can also account for sym¬ 
metric outer velocity warps, and can correct for mi¬ 
nor distortions due to seeing. Unlike other methods, 
DiskFit fits a specific physically motivated model rather 
tha n parametrizing concentric rings of the velocity field 


tua n parametrizing concentric rings ot tue velocity Ueld 
(see S£ekkens_^_SelIwood 20071 : ISellwood fc Sanched I 2 OI 0 I : 
iKuzio de Narav et al.ll2012l i . DiskFit also returns realistic 
estimates of the uncertainties on the best fitting parameters 


^ http://www.physics.rutgers.edu/~spekkens/diskfit/ 


using bootstrap realizations of the best fitting model. Disk- 
Fit is therefore ideally suited for direct searches for bar-like 
flows in nearby galaxies. 

Large resolved spectroscopic surveys have provided a 
method for statistically sampling the galaxy population, al¬ 
lowing for a greater understanding of their formation and 
evolution. The large number of objects that are studied al¬ 
lows for a meaningful statistical analysis over a wide range 
of galaxy types and environmental conditions. This has led 
recently to the advent of several large kinematic studies 
using Integral Field Units (IFUs). Early IFU surveys typ¬ 
ically targe ted 10 — 100 near by galaxies, and include 
SAU RON dde Zeeuw et al.l I2OO2II. VENG A dBlanc et al.l 
I2OI3II. PINGS dRosales-Ortega et al.ll2010ll . and ATLAS3D 
( Cappel lari et al.| I2OII ). Ut )coming large surveys include 


MaNGA dBundv et al 


__ _ ., . J| 2015fl as w ell as the ongoing SAMI 

dCroom et al.l ^12l: Br yant et al']l2015l l and CALIFA sur¬ 


veys (Sanchez et al. IM). 

The CALIFA (Calar Alto Legacy Integral Field Spec¬ 
troscopy Area) Survey utilizes the PPaK mode of the PMAS 
IFU on the 3.5 m telescope at the Calar Alto observatory 
dRoth et al.l l2005l : iKelz et ^ I 2 OO 6 II and will produce spa¬ 
tially resolved spectral information for nearly 600 galaxies. 
The CALIFA mother sample consists of 939 possible target 
galaxies and was selected fro m the SPSS Data Rel ease 7 
(DR7) photometric catalogue dAbazaiian et al]l2009ll . Tar¬ 
gets are drawn according to visibility from this mother sam¬ 
ple. The fully calibrated, first data release (DRl) DRl of 
100 CALIFA targets occurred on 1 November 2012 and is 
the source of the data used in this paper. 

Although CALIFA will produce a much smaller sample 
than either MaNGA or SAMI, the number of resolution el¬ 
ements that covers each target galaxy is much greater. The 
CALIFA survey is therefore unique in its combination of res¬ 
olution and sample size when compared both to past (e.g. 
PINGS, VENGA) and future (e.g. SAMI, MaNGA) spec¬ 
troscopic surveys. This resolution is particularly useful in 
searching for non-circular flows, where several resolution el¬ 
ements across both the major and minor axes are required. 

Several studies of the kinematic properties of barred 
galaxies have now been carried out by the CALIFA col¬ 
laboration. The ionized gas kinematics for a large sam¬ 
ple of CALIFA galaxies with a wide range of mo rpho¬ 
logical types were examine d by iKehrig et al.l d2012ll and 
ICarcfa- Lorenzo et al.l d2014l L The pa ttern speeds of barre d 
CALIFA galaxies were examined by lAguerri et al.l d2015ll , 
and the kinematic align ment of barred and unbar r ed sys - 
tems were compared by iBarrera-Ballesteros et ah 1 (I20l4) . 
However, no attempt to separate non-circular flows from the 
underlying disk rotation has yet been made. 

This paper carries out the first direct search for non¬ 
circular motions in the H a velocity fields of an intermediate- 
inclination subsample of CALIFA DRl galaxies using Disk- 
Fit. The kinematic bar search performed in this paper pro¬ 
vides a direct probe of the dynamical impact of the bar on 
the disk - and thus insight on its role in driving galaxy evo¬ 
lution. In addition, a systematic comparison between kine¬ 
matic and photometric bar classifications is carried out, af¬ 
fording both a preliminary look at the relationship between 
these two schemes as well as a search for systems where the 
photometric classification belies a different kinematic struc¬ 
ture. This paper therefore provides a first look at the re- 
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lationship between photometric and kinematic indicators of 
bar-like flows in disk galaxies. 


2 SAMPLE SELECTION 

We attempt to model all intermediate in clination gas rich 
CALI FA DRl galaxies. As des cribed in iHusemann et all 
(l2013ll and IWaldier et ^ (12014) . DRl is a statistically rep¬ 
resentative subset of the full CALIFA sample in stellar mass, 
colour, and morphology. 

The morphological classification and barredness of all 
CALIFA galaxies wer e determined by-eye by five members 
of the collaboration dWalcher et al.lT20l4 l. In the sample 
of 100 CALIFA DRl galaxies, 30% are classified as barred 
(B), 21% as intermediately barred (AB), and 49% as un¬ 
barred (A). These values are consistent with the optical bar 
fractions in nearby galaxies discussed in !JT] Figure [T] shows 
SDSS 2' X 2' igr composite images of galaxies in the CALIFA 
mother sample illustrating the different bar classifications. 
We adopt these morphological and bar classifications for our 
analysis. The photometric inclinations, ip, for each galaxy 
were computed from ellipses enclosing 50-90% of the total 
galaxy flux in its r-band SDSS image. We attempt to model 
only those DRl galaxies with 40° < ip < 75°. At lower incli¬ 
nations, the line-of-sight component of the rotation velocity 
becomes harder to disentangle from the turbulent motions 
in the disk, an d the inclination is co rrespondingly harder to 
constrain fe.g. iBershadv et al.ll201^ . It is possible to mea¬ 
sure rotation curves at higher inclinations, but searching for 
non-circular flows becomes difficult because there is limited 
minor axis information. 57/100 of the CALIFA DRl galaxies 
meet this photometric inclination selection criterion. 

The H a velocity fields of these 57 galaxies were exam¬ 
ined to determine their suitability for modelling. The Ho? 
emission line maps were extracted from the gas-pure cubes 
by fitting a single Gaussian function using the systemic ve¬ 
locity from the anal ysis of the stell ar population as an initial 
guess using FIT3D (ISanchej|2006l '). The uncertainties in the 
velocity and velocity dispersion were determined from the in¬ 
verse of the co-variance map. The velocity fields used could 
be affected by the presence of dust by both hindering the 
detection of the emission lines as well as potentially perturb - 
ing the velocity distribution dReean. Sheth fc Vogell[l9^ . 
The obscuration of Ha emission could affect our search, if 
there is a significant lack of Ha along the major or minor 
axis of the galaxy, a search for non-circular flows would not 
be possible (see >14.11) . Typically the presence of strong dust 
lanes is associated with spir al arms or the rotational front 
of bars d Athanassoulal Il992l f and is unlikely to produce or 
mimic the coherent motion for which we are searching, and 
therefore would not significantly impact our results. 

We apply spatial masks created from SDSS r-band 
masks reprojected to the scale, orientation and pixel size 
of the CALIFA cubes with foreground stars and artefacts 
removed. All those galaxies where less than 5% of the pixels 
contained H a were eliminated from the sample, and we at¬ 
tempted to model galaxies where 5 — 10% of them contained 
H a only if a velocity gradient was discernible in the velocity 
field. This criterion eliminated nearly all remaining ellipti¬ 
cal galaxies, resulting in a sample of 49/100 CALIFA DRl 
galaxies. We attempted to model each of these 49 galaxies. 


The basic properties of these 49 galaxies are given in 
Table [U and Figure [2] shows their distribution of morpho¬ 
logical types and bar classification. In Figure [51 the dotted 
histogram shows the distribution of all 100 CALIFA DRl 
galaxies and the blue histogram shows the 49 galaxies that 
we attempted to model. The modelled galaxies span a range 
of morphological types as well as barredness, in similar pro¬ 
portions to the full DRl sample. The red histograms in 
Figure[5]show the galaxies for which valid kinematic models 
were obtained; we describe our process for selecting valid 
models in >14.11 


3 KINEMATIC MODELS 

We use the publicly available DiskFit code to carry out all 
of our kinematic modelling. DiskFit returns the galaxy cen¬ 
tre, kinematic inclination {ik), systemic velocity (Vi,;/s), and 
position angle (PA) of the kinematic major axis. We con¬ 
sider three types of kinematic models within DiskFit. The 
simplest model includes rotation only, and is given by: 

Fmodei (f^) — Fsys “t“ Ft (^) siu ik COS 0 (1) 

where Vays is the systemic velocity, Vt is the rotation ve¬ 
locity, ik is the kinematic disk inclination, and 9 is the az¬ 
imuthal angle from the major axis in the plane of the disk. 
Note that DiskFit assumes that the disk is flat, and returns 
a single value for the position angle, inclination, systemic 
velocity, and disk centre. 

The bisymmetric model introduces non-circular flows 
produced by a bar-like, m = 2 perturbation to the potential 
and is given by: 

^modellFCj — Vsys -f Siu ^Vf(R)cOS0 (2) 

— V2,t{R) COS {29b) COS 9 — V2,r{R) sin {29b) sin0] 

where V2,t and Vi.r are the tangential and radial components 
of the non-circular flow respectively and the major axis of 
the bar is at an angle 9b to the projected major axis in 
the plane of the disk. Equation [2] describes first-order gas 
flows in a barred galaxy. Although higher-order harmonic 
components are required to fully describe the kinematics 
of strongly barred galaxies, the m = 2 component always 
dominates and is generally sufficient to model the physical 
properties of the flow. 

After an initial search for bar-like flows across the entire 
disk extent in all galaxies, an outer search radius (Rmaa;) was 
imposed, beyond which all non-circular flows are forced to 
zero. For all galaxies with bar classification A or AB, this 
radius is one half that of the disk. For those galaxies with 
bar classification B, Umax is the larger of either one half the 
galaxy radius or the estimated bar radius plus 20% to ensure 
that the entire bar is included. 

We also applied a model with radial flows that is also 
available in DiskFit and assumes an m = 0 perturbation in 
the disk plane: 

Vmodei{R) = VsyB{R) + sini [Vt{R) cos9 + Vr{R) sin6>] (3) 

where K- is the radial velocity. Since quiescent disk galaxies 
are unlikely to exhibit large radial flows, we focus our anal¬ 
ysis on the rotation only and bisymmetric DiskFit models. 

DiskFit can also fit for a warp in the outer disk. The 
disk is assumed to be flat out to some radius after which 
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Table 1. Literature values for 49/100 CALIFA DRl galaxies that were modelled with DiskFit. 


ID 

(1) 

Name 

(2) 

Morph 

(3) 

Bar 

(4) 

a 

(hh:mm:ss.s) 

(5) 

5 

(dd:mm:ss.s) 

(6) 

ip 

(deg) 

(7) 

(deg) 

(8) 

Lsys 

(km s~^) 

(9) 

Ref. 

(10) 

Mod. 

(11) 

3 

NGC7819 

Sc 

A 

00:04:24.5 

31:28:19.2 

59 

65 

4958±6 

1,1 

Y 

7 

UGC36 

Sab 

AB 

00:05:13.9 

06:46:19.3 

53 

20 

6291±13 

1,2 

Y 

10 

NGC36 

Sb 

B 

00:11:22.3 

06:23:21.7 

49 

170 

6030±4 

1,3 

Y 

14 

UGC312 

Sd 

B 

00:31:23.9 

08:28:00.2 

72 

5 

4364±3 

1,2 

N (DF) 

42 

NGC477 

Sbc 

AB 

01:21:20.5 

40:29:17.3 

49 

120 

5876±8 

1,3 

Y 

43 

IC1683 

Sb 

AB 

01:22:38.9 

34:26:13.7 

54 

160 

4869±20 

1,21 

N (DF) 

73 

NGC776 

Sb 

B 

01:59:54.5 

23:38:39.4 

46 

135 

4921±6 

1,4 

N(I) 

100 

NGC1056 

Sa 

A 

02:42:48.3 

28:34:27.0 

56 

155 

1545±10 

1,4 

Y 

146 

UGC3253 

Sb 

B 

05:19:41.9 

84:03:09.4 

51 

65 

4130±7 

1,4 

Y 

156 

NGC2449 

Sab 

AB 

07:47:20.3 

26:55:48.7 

61 

138 

4892±3 

5,6 

Y 

277 

NGC2916 

Sbc 

A 

09:34:57.6 

21:42:18.9 

54 

18 

3730±7 

5,4 

Y 

306 

UGC5358 

Sd 

B 

09:58:47.1 

11:23:19.3 

74 

75 

2914±1 

5,7 

N (DF) 

307 

UGC5359 

Sb 

B 

09:58:51.6 

19:12:53.9 

70 

94 

8472±9 

5,8 

Y 

341 

UGC5771 

E6 

A 

10:37:19.3 

43:35:15.3 

45 

52 

7403±52 

5,9 

Y 

479 

NGG4003 

SOa 

B 

11:57:59.0 

23:07:29.6 

67 

167 

6509±32 

5,4 

N (DD) 

486 

UGC7012 

Scd 

AB 

12:02:03.1 

29:50:52.7 

58 

10 

3081±1 

1,10 

Y 

515 

NGC4185 

Sbc 

AB 

12:13:22.2 

28:30:39.5 

50 

164 

3904±5 

5,11 

Y 

518 

NGC4210 

Sb 

B 

12:15:15.8 

65:59:07.2 

43 

100 

2732±7 

5,3 

Y 

528 

IC776 

Sdm 

A 

12:19:03.1 

08:51:22.2 

56 

98 

2468±1 

4,12 

Y 

548 

NGG4470 

Sc 

A 

12:29:37.8 

07:49:27.1 

49 

1 

2341±1 

5,12 

Y 

608 

NGG5000 

Sbc 

B 

13:09:47.5 

28:54:25.0 

53 

6 

5608±4 

5,4 

N(I) 

657 

UGC8733 

Sdm 

B 

13:48:39.0 

43:24:44.8 

61 

140 

2338±6 

1,4 

Y 

665 

UGC8781 

Sb 

B 

13:52:22.7 

21:32:21.7 

59 

175 

7592±15 

1,13 

Y 

676 

NGG5378 

Sb 

B 

13:56:51.0 

37:47:50.1 

51 

78 

3042±25 

5,14 

Y 

680 

NGC5394 

Sbc 

B 

13:58:33.2 

37:27:13.1 

41 

60 

3448±2 

5,6 

N (DF) 

764 

NGC5720 

Sbc 

B 

14:38:33.3 

50:48:54.9 

49 

52 

7790±6 

5,9 

Y 

769 

UGC9476 

Sbc 

A 

14:41:32.0 

44:30:46.0 

51 

132 

3262±8 

5,3 

Y 

823 

NGG6063 

Sbc 

A 

16:07:13.0 

07:58:44.4 

53 

155 

2848±5 

5,15 

Y 

824 

IC1199 

Sb 

AB 

16:10:34.3 

10:02:25.3 

67 

159 

4731±44 

5,9 

Y 

826 

NGC6081 

SOa 

A 

16:12:56.9 

09:52:01.6 

64 

129 

5176±31 

5,9 

Y 

833 

NGC6154 

Sab 

B 

16:25:30.5 

49:50:24.9 

50 

166 

6015±40 

5,4 

Y 

846 

UGC10695 

E5 

A 

17:05:05.6 

43:02:35.4 

47 

117 

8328±34 

5,22 

N (DD) 

850 

NGG6314 

Sab 

A 

17:12:38.7 

23:16:12.3 

60 

174 

6633±4 

5,16 

Y 

852 

UGC10796 

Scd 

AB 

17:16:47.7 

61:55:12.4 

62 

62 

3079±11 

5,5 

N (DF) 

854 

UGC10811 

Sb 

B 

17:18:43.7 

58:08:06.4 

67 

92 

8746±26 

5,17 

Y 

856 

IC1256 

Sb 

AB 

17:23:47.3 

26:29:11.5 

54 

90 

4730±10 

1,4 

Y 

858 

UGC10905 

SOa 

A 

17:34:06.4 

25:20:38.3 

59 

0 

7843±34 

1,18 

Y 

863 

NGC6497 

Sab 

B 

17:51:18.0 

59:28:15.2 

50 

116 

6162±64 

5,9 

Y 

865 

UGC11228 

SO 

B 

18:24:46.3 

41:29:33.8 

56 

165 

5771±23 

1,23 

N (DF) 

866 

UGC11262 

Sc 

A 

18:30:35.7 

42:41:33.7 

69 

48 

5606±36 

4,18 

Y 

867 

NGG6762 

Sab 

A 

19:05:37.1 

63:56:02.8 

62 

115 

2923±47 

1,9 

Y 

874 

NGG7025 

SOa 

A 

21:07:47.3 

16:20:09.2 

44 

50 

4968±5 

1,19 

Y 

877 

UGC11717 

Sab 

A 

21:18:35.4 

19:43:07.4 

63 

35 

6303±36 

1,9 

Y 

887 

NGC7321 

Sbc 

B 

22:36:28.0 

21:37:18.4 

46 

25 

7145±5 

1,19 

Y 

890 

UGC12185 

Sb 

B 

22:47:25.1 

31:22:24.7 

63 

150 

6649±10 

1,4 

Y 

896 

NGG7466 

Sbc 

A 

23:02:03.5 

27:03:09.3 

59 

25 

7508±3 

1,20 

Y 

901 

NGC7549 

Sbc 

B 

23:15:17.3 

19:02:30.4 

41 

110 

4736±3 

1,24 

N (DF) 

904 

NGC7591 

Sbc 

B 

23:18:16.3 

06:35:08.9 

54 

170 

4956±4 

1,4 

Y 

935 

UGC12864 

Sc 

B 

23:57:23.9 

30:59:31.5 

74 

80 

4683±7 

1,4 

N (DF) 


Notes. Col. (1): CALIFA ID. Col. (2): Galaxy name. Col. (3): Morphological type, as classified by CALIFA. Col. 
(4): CALIFA bar classification. Col. (5): Photometric right ascension as provided by NED. Col. (6): Photometric 
declination as provided by NED. Col. (7): Photometric inclination. Col. (8): Photometric disk PA. Col. (9): 
Systemic heliocentric velocity. Col. (10): Reference key for (pp, and Vsys. Col. (11): Y: Successfully modelled. N: 
Unable to model, not included in final sample (DD: distorted inner disk. I: inclination out of range. DF: not in 

dynamical equilibrium). _ _ _ _ _ _ 
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The incidence of bar-like kinematic flows in CALIFA galaxies 5 



Figure 1. SDSS 2' x 2' igr composite images of galaxies in the CALIFA mother sample illustrating the different bar classifications, (a) 
A: NGC 1056. (b) AB: NGC 4185. (c) B: NGC 7321. 
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Figure 2. Distribution of morphological type (left) and photometric bar strength (right) for the galaxies modelled with DiskFit. In 
both panels, the dotted line shows the distribution of all 100 CALIFA DRl galaxies, the blue histogram shows the distribution for the 
subset of 49 intermediate-inclination, gas-rich galaxies that we attempted to model with DiskFit (see and the red histogram shows 
the distribution of our final sample of 37 galaxies for which valid kinematic models were obtained (see S- The modelled galaxies span 
a range of morphological type and barredness. 


the ellipticity and PA vary quadratically with increasing ra¬ 
dius. There are degeneracies between fitting for a warp and 
bisymmetric flows since both can cause variations in the el¬ 
lipticity and PA of the flow pattern, and therefore both of 
these components cannot be fit simultaneously. We apply a 
warp model in only a few cases where non-circular flows are 
evident but not well parametrized by other models. 

Multiple DiskFit models were attempted for the 49/100 
DRl galaxies selected using the criteria in 32] DiskFit was 
applied to each galaxy automatically using a series of scripts 
to drive the publicly available executable. We adopt ip de¬ 
scribed In ^as an initial guess for it and use the literature 
Vaya values from Table [T] as an Initial guess for that parame¬ 
ter. The PA and Rmax were visually estimated from contour 
plots of the input velocity field, and the centre position was 
initially chosen as the central pixel of the PMAS IFU. Al¬ 
though the atmo spheric seeing in the CALIFA observations 
is typically ~ 1" dSanchez et al.ll2012l l the spatial resolution 
of the velocity field is ~ 3.5”. The ring radii at which the 


velocity components were sampled by DiskFit were spaced 
by 3”. 


Initially rotation only (Eq. [!}, bisymmetric (Eq. [2]), and 
radial (Eq. [2| models were applied to the velocity fields for 
each of the 49/100 galaxies selected as described in 32] We 
generate 100 bootstrap realizations of each velocity field 
to deter mine uncertainties on each m odel parameter as de¬ 
scribed in ISellwood fc Sanchez! d2010l L The parameter AISM 
is added in quadrature to the uncertainties in the emission 
line centroids during the fit w hich allows for turbulence i n 
the interstellar medium (ISM) (ISoekkens fc SellwoodllioOTl l . 
A value of AISM= 5kms“^ was used for all models. For all 
successfully modelled galaxies (see Q, uncertainties on all 
parameters were derived from 1000 bootstrap realizations of 
the best-fitting model. 
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4 RESULTS 

In the following sections we detail the process by which an 
optimal DiskFit model was selected for each of the modelled 
galaxies (ED), provide an in-depth case study for two rep¬ 
resentative systems (NGC 1056 and NGC 7321, >14.211 . and 
present the results from the DiskFit models for the sample 
as a whole M- 


4.1 Model Selection 


Each of the DiskFit models described in i}^ was applied to 
the 49/100 DRl galaxies that satisfy the criteria described 
in ^ We then proceeded to determine which model pro¬ 
vided the best characterization of the data. The overall fit 
quality was assessed by examining the model, (data - model) 
residuals, disk geometries, and rotation curves for the rota¬ 
tion only, bisymmetric, and radial models. Since the radial 
model was used primarily for diagnostic purposes, we fo¬ 
cused on determining whether the rotation only or bisym¬ 
metric model provided the best description of each system. 
We illustrate the details of this process for two representa¬ 
tive sample galaxies in >14.21 

In addition to examining the DiskFit outputs described 
above, a chi-square test was used to determine the statis¬ 
tical significance of any non-circular flows returned by the 
bisymmetric and radial model fits. The statistic for the 
null hypothesis that the non-circular flows are consistent 
with zero is given by: 


N 

i=l 


(xi - fii Y 


E' 


(4) 


where Xi is the velocity (14,t or 14,r for the bisymmetric 
model), CTi is its uncertainty, N is the number of radial bins 
in which the flow is measured, and = 0. We deem that 
values of Xnc > Scr for the bisymmetric model, where a 
is the standard deviation of the distribution with N de¬ 
grees of freedom, signal statistically significant flows, values 
3'^ < X%c < represent marginally significant flows, and 
insignificant flows have Xnc < 3cr. The chi-square test was 
found to be useful in determining the best physical model 
but was not always accurate. Specifically, the significance 
threshold of 5cr was too stringent in four cases (NGC 477, 
UGC 3253, NGC 5720, and UGC 12185) where the uncer¬ 
tainties are large for the bisymmetric model. Of these, in 
two cases (NGC 5720 and UGC 12185) Xnc > if 
centre position is fixed, and all four had Xnc > 5a if the 
radial model was examined. The chi-square test is therefore 
a good indicator of significant non-circular flows, but it was 
still necessary to examine by eye the residuals, disk geome¬ 
try and kinematic components for each model to determine 
which was optimal. 

The optimal physical model for each galaxy fell into 
one of two possible categories: Rotation Only (25/37 or 
67.6itg|% galaxies), and Non-Circular Flows (12/37 or 
32.4lg g% galaxies). Note that the uncertainties on derived 
fractions correspond to the Icr c onfidence i n terva ls obtained 
using the Bayesian approach of ICamero 3 lloTJ). The Ro¬ 
tation Only category includes all galaxies where either a 
rotation only model was deemed to be optimal {17125 or 
68 . 0 /j/^ 4 % galaxies), or where a search for bisymmetric 
flows was not possible (8/25 or 32.0/7'^7^% galaxies). This 


Can’t Tell subcategory contains galaxies with velocity fields 
lacking significant H a emission or containing masked pixels 
near the galaxy centre. The majority of the galaxies that fell 
into the Non-Circular Flows category have significant non¬ 
circular flows in the bisymmetric model that are consistent 
with bar-like flows. We return to this issue in >J3 

Out of the 49/100 galaxies that we attempted to model, 
10 galaxies were eliminated from further consideration be¬ 
cause they did not meet the assumptions inherent in the 
DiskFit models. Two galaxies (NGG 4003 and UGG 10695) 
had model residuals that suggested they have strongly dis¬ 
torted inner disks, and eight were determined to be dis¬ 
turbed and unlikely to be in dynamical equilibrium as im¬ 
plicitly assumed by any model that fits for a rotating disk. 
Although fixing parameters duri ng kinematic model fitting 
is common in the literature (e.g. Ide Blok et al.l[2003 ') it was 
decided that if DiskFit could not find the inclination of 
the galaxy, even in the simplest rotation only model, that 
it was unlikely that the galaxy could be well-modelled as 
an equilibrium fiat disk. There are therefore 39/100 GAL- 
IFA DRl galaxies for which valid models were determined. 
Table [T] indicates which of the 49/100 galaxies were success¬ 
fully modelled and those that were eliminated for these rea¬ 
sons. For consistency with the photometric selection criteria 
in >J 3 we eliminated an additional two systems for which 
ik < 40°. Our final sample therefore consists of 37/100 GAL- 
IFA DRl galaxies, whose distribution in morphological type 
and barredness is shown by the red histograms in Figure (2] 
The results of the best fitting models for the final sam¬ 
ple galaxies are shown in Table [2] The velocity field, veloc¬ 
ity uncertainties, best fitting model and residuals for each 
galaxy can be found online in Appendix A. The H a velocity 
fields for the 12 galaxies where fits were attempted but then 
rejected can be found online in Appendix B. Detailed notes 
on all of the 49/100 galaxies are online in Appendix C. 


4.2 Model Selection for Representative Galaxies 

In this section, we illustrate the optimal DiskFit model se¬ 
lection process outlined in >14. II for the representative galax¬ 
ies NGC 1056, a galaxy that we classify as having rotation 
only flows, and NGC 7321, which we classify as having non¬ 
circular flows. For both galaxies, the best fitting disk pa¬ 
rameters for the rotation only and bisymmetric models are 
given in Tabled 


A.2.1 NGC 1056 

NGC 1056 is at a distance of 18.7 Mpc dTheureau et al.l 
I 2 OO 7 II . It is classified as having m orphological type Sa with 
photometric bar classification A JWalcher et al.l[2014l l and 
is a typical example of a galaxy whose Ha velocity field 
is best described by a rotation only model. Panel (a) of 
Figure [ 1 ] shows the SDSS 2' x 2' igr composite image of this 
galaxy. Figure Oshows the H a velocity field and correspond¬ 
ing uncertainties, as well as the best-fitting rotation only 
and bisymmetric models and their residuals. There is Ha 
emission throughout the disk and the velocity field is well 
determined, with uncertainties Avobs < 7 kms“^ through¬ 
out most of the inner disk. The residuals for the rotation 
only model show no apparent coherent patterns, and the 
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Table 2. DiskFit minimization results for final sample of 37 galaxies. 


ID 

(1) 

Name 

(2) 

Model 

(3) 

Sig 

A) 

(4) 

X 

(arcsec) 

(5) 

y 

(arcsec) 

(6) 

(deg) 

(7) 

'^k 

(deg) 

(8) 

Vsys 

(km s“^) 

(9) 

<t>b 

(deg) 

(10) 

(deg) 

(11) 

3 

NGC7819 

m2 

5.5 

l.liO.l 

-1.2±0.1 

80±2 

45±4 

4955±1 

110±6 

18,95 

7 

UGC36 

R 

0.2 

l.OiO.l 

-1.8±0.1 

20±1 

65±2 

6303±2 



10 

NGG36 

W 

0.8 

0.9±0.3 

-1.8±0.6 

200±1 

51±2 

5992±5 



42 

NGG477 

m2 

1.8 

5.1±0.4 

lO.OiO.3 

311±2 

38±5 

5875±3 

19±38 

326,244 

100 

NGG1056 

R 

1.3 

1.4±0.2 

-1.9±0.2 

161±1 

58±2 

1559±1 



146 

UGC3253 

m2 

1.8 

0.9±0.2 

-1.7±0.3 

84±1 

52±3 

4119±1 

59±11 

64 

156 

NGG2449 

CT 

1.7 

1.4±0.2 

-1.3±0.2 

136±1 

58±1 

4903±2 



277 

NGG2916 

CT 

1.9 

0.4±0.2 

0.2±0.7 

16±1 

51±1 

3698±5 



307 

UGC5359 

R 

1.6 

1.2±0.1 

-l.OiO.l 

273±1 

69±1 

8465±1 



341 

UGC5771 

CT 

0.1 

1.0±0.2 

-1.1±0.2 

54±1 

34±8 

7417±4 



486 

UGC7012 

R 

0.5 

1.5±0.2 

-0.3±0.6 

6.±1 

52±4 

3088±3 



515 

NGG4185 

R 

0.0 

1.0±0.2 

-1.4±0.2 

168±1 

48±1 

3873±1 



518 

NGG4210 

CT 

0.0 

0.9±0.7 

-0.9±0.3 

97±1 

41±1 

2712±3 



528 

IC776 

R 

0.8 

1.9±0.5 

0.5±0.6 

281±1 

46±4 

2470±1 



548 

NGG4470 

R 

0.0 

0.5±0.7 

0.2±0.4 

171±1 

47±2 

2344±1 



657 

UGC8733 

CT 

1.5 

8.1±1.3 

0.7±1.2 

218±3 

63±5 

2335±5 



665 

UGC8781 

CT 

4.1 

1.2±0.2 

-1.3±0.5 

160±1 

58±2 

7571±4 



676 

NGG5378 

CT 

0.1 

-0.2±0.3 

-0.9±0.3 

84±2 

39±5 

2960±2 



764 

NGG5720 

m2 

8.3 

1.2* 

-1.0* 

309±1 

51±1 

7785±1 

47±4 

278 

769 

UGC9476 

R,ad 

0.0 

2.1±0.4 

-1.2±0.3 

117±3 

43±2 

3247±2 



823 

NGG6063 

R 

0.1 

1.1±0.2 

-0.9±0.1 

152±1 

57±1 

2841±1 



824 

IC1199 

R 

0.7 

0.8±0.1 

-1.6±0.2 

158±1 

61±2 

4708±2 



826 

NGG6081 

R 

0.1 

1.7±0.4 

-1.5±0.3 

129±1 

68±3 

5050±5 



833 

NGG6154 

CT 

0.1 

0.4±0.3 

O.liO.3 

202±3 

70±6 

5983±2 



850 

NGG6314 

R 

1.9 

1.6±0.1 

-0.5±0.2 

176±1 

62±1 

6614±4 



854 

UGC10811 

m2 

6.9 

0.7±0.2 

-1.6±0.1 

90±1 

69±1 

8739±3 

64±5 

125 

856 

IC1256 

R 

0.0 

1.4±0.2 

-0.8±0.1 

270±1 

52±1 

4717±1 



858 

UGC10905 

R 

0.1 

l.liO.l 

-0.5±0.2 

175±1 

60±1 

7750±3 



863 

NGG6497 

m2 

6.6 

1.2±0.2 

-1.0±0.2 

114±1 

57±1 

6053±1 

60±8 

157 

866 

UGC11262 

R 

0.0 

0.7±0.1 

-0.9±0.1 

54±1 

68±1 

5546±1 



867 

NGG6762 

R 

5.2 

1.7±0.2 

-1.2±0.2 

299±2 

64±4 

2939±3 



874 

NGG7025 

m2 

14.0 

l.OiO.l 

-l.OiO.l 

40±1 

62±1 

4925±2 

68±2 

89,29 

877 

UGC11717 

R 

0.4 

0.7±0.4 

-1.2±0.4 

224±1 

59±2 

6272±8 



887 

NGG7321 

m2 

9.7 

1.6±0.2 

-1.4±0.3 

12±1 

46±2 

7123±3 

47±6 

48 

890 

UGC12185 

m2 

7.1 

0.6* 

1 

o 

337±1 

59±1 

6586±1 

44±5 

310 

896 

NGG7466 

R 

1.2 

0.6±0.2 

-1.2±0.3 

25±1 

62±2 

7489±3 



904 

NGG7591 

m2 

13.0 

0.6±0.1 

-l.liO.l 

145±1 

60±2 

4929±1 

64±4 

190 


Notes. Col. (1): CALIFA ID. Col. (2): Galaxy name. Col. (3): Optimal DiskFit model: R (rotation only); CT 
(rotation only, can’t tell); m2 (bisymmetric); Rad (radial); W (warp). Col. (4): Maximum a value from chi-square 
test on significance of non-circular flows (Vi.t or V 2 ,r) in bisymmetric model. Col. (5): Right ascension of disk 
centre, relative to photometric center from Table [J * = centre was fixed. Col. (6): Declination of disk center, 
relative to photometric center from Table [T] * = centre was fixed. Col. (7): Disk PA. Col. (8): Disk inclination. 
Col. (9): Disk systemic velocity. Col. (10): Bisymmetric distortion PA in the disk plane. Col. (11): Bisymmetric 
distortion PA in the sky plane. If no photometric bar is evident in the galaxy image, possible values for either the 
major or minor axis are listed. 


bisymmetric model does not significantly reduce the resid¬ 
uals despite the additional free parameters in the fit and 
lower reduced xlt relative to the rotation only fit. As shown 
in Table [S] there is good agreement between the derived 
kinematic parameters for both models. Note that the rel¬ 
atively high values of Xa for both models suggested that 
AISM= 5kms“^ is too low for this galaxy. 

The rotation only model was chosen for NGC 1056 be¬ 
cause there is no evidence for non-circular flows in the model 
residuals. A chi-square test run on the components of the 
bisymmetric flow (Eq. supports this assessment, which 
returned xl/c < 3 (j and therefore no evidence against the 
null hypothesis of a purely rotating disk. 

NGC 1056 is one of few sample galaxies where a pub¬ 


lished rotation curve derived from an independent H a veloc¬ 
ity field is available for comparison. Figure [4] shows the Disk- 
Fit rotation curve for the rotat ion only model (blue curve) 
compa red to that published bv lEpinat, Amram fc MarcelinI 
ll2008l) from the GHASP survey (solid red and green 
curves for the approaching and receding sides respec¬ 
tively). The published rotation curve has an amplitude 
of Vt ~ 160 kms“^ which is significantly greater than 
that found with DiskFit of Vt ~ 130 kms“^. This dif¬ 
ference arises primarily due to the difference in kine- 
matic inclination found by DiskF it (ik = 58 ± 2°) and by 
lEpinat. Amram fc MarcelinI ll20 08l) (ft. = 41 ± 1 0 °V I'O" 
tation curve from Epinat. Amram fc MarcelinI (l2008l ) was 
therefore reprojected using the DiskFit inclination, and 
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Table 3. DiskFit minimization results for representative galaxies NGC 1056 and NGC 7321. 


ID 

(1) 

Galaxy 

(2) 

X 

(arcsec) 

(3) 

y 

(arcsec) 

(4) 

(deg) 

(5) 

(deg) 

(6) 

Vsy s 

(km s~^) 

(7) 

N 

(8) 

(9) 

100 

NGC1056 
Rotation Only 

0.1±0.2 

-1.9±0.2 

161±1 

58±2 

1559±1 

2781 

3.8 


Bisymmetric 

0.1±0.2 

-1.9±0.2 

161±1 

59±1 

1559±1 

2776 

3.5 

887 

NGC7321 
Rotation Only 

O.liO.l 

-1.3±0.3 

10±1 

45±3 

7122±3 

2837 

1.5 


Bisymmetric 

0.1±0.2 

-1.4±0.3 

12±1 

46±2 

7123±3 

2822 

0.9 

Notes. 

Col. (1): CALIFA ID. 

Col. (2): Galaxy 

name and mo 

dels. Col. (3): 

Right ascension of disk centre, relative 


to photometric center from Table ITI Col. (4): Declination of disk center, relative to photometric center from Table 
m Col. (5): Disk PA. Col. (6): Disk inclination. Col. (7): Disk systemic velocity. Col. (8): Number of degrees of 
freedom. Col. (9): Model reduced chi-square. The optimal model for each galaxy is shown in bold. 


0 


r (kpc) 
2 



r (orcsec) 


Figure 4. Comparison between DiskFit rotation curve 
for NGC 1056 and that from the GHASP survey 
llEpinat. Amram &: MarcelinI l20Q8h . The blue curve is the 
rotation curve from the rotation only model shown in 
Figur e!^ The red and green s olid curves are velocities 
from lEpinat, Amram &: MarcelinI bOOSl l for the approach¬ 
ing and receding sides respectively. The dashed red and 
green curves are the approaching and receding velocities from 
lEpinat, Amram &: MarcelinI bonSil rescaled from their kinematic 
inclination of jj. = 41° to the DiskFit-derived ij. = 58°. 


shown as the dashed red and green curves in Figure S) 
There is reasonable agreement b etween the results from 
lEpinat. Amram fc MarcelinI (l2008l i and DiskFit’s rotation 
only model when the same disk inclination is adopted in 
both models in the inner 5” (receding) and beyond 15” (ap¬ 
proaching), considering that Epinat et al. modelled differ¬ 
ent data using a differ ent algorithm, and considered only 
half the disk at a time. lEpinat. Amram fc MarcelinI (l2008l l 
also used adaptive binning techniques to generate the ve¬ 
locity field from which their rotation curve was derived, 
which may explain the discrepancies between their results 
and ours for 5” < r < 15”. In general, the disk pa¬ 
rameters found by DiskFit (Table Et agree well with those 
from the literature (Table [H Barrera-Ballesteros et al.ll20l4 

iGarcfa-Lorenzo et al. lE^. 


4.2.2 NGC 7321 

NGC 7321 is at a distance of 97.5 Mpc (|Theureau et al.l 
l2007h . It is classified as having mo rphological type Sbc with 
photometric bar classification B (jWalcher et aklEoilh and 
is a typical example of a galaxy whose Ha velocity field 
is best described by a bisymmetric model. Panel (c) of 
Figure [T] shows the SDSS 2' x 2' igr composite image of 
this galaxy. Figure [5] shows the Ha velocity field and cor¬ 
responding uncertainties, as well as the best-fitting rotation 
only and bisymmetric models and their residuals. Uncertain¬ 
ties in the disk are all relatively small with Avobs < 7 kms“^ 
across most of the disk, with only small regions of higher un¬ 
certainties, up to Avobs ~ 25 kms“^. Note that while there 
is a large area masked in the top left of the velocity field 
due to a foreground star, it does not strongly affect the fit 
because of its location with respect to the major and minor 
axis. 

Although the rotation only model (Figure Eh) fits the 
data reasonably well, the velocity field exhibits ‘S-shaped’ 
isovelocity contours along the minor axis that (by definition) 
are not replicated in this model. As a result, a coherent fea¬ 
ture (residuals negative for 0 < x < 10” east of the galaxy 
centre, and positive residuals 0 < x < 10” to the west) is 
present near the centre of the residual map. The bisym¬ 
metric model and residuals are shown in the bottom row 
of Figure m The model appears to be a good fit, similar 
to the rotation only model, except it now accounts for the 
‘S-shape’ feature along the minor axis. A significant differ¬ 
ence is noted when comparing the residuals of the rotation 
only and the bisymmetric models: the coherent rotation only 
model residual pattern is not present in the bisymmetric 
model residuals. The kinematic components for the bisym¬ 
metric model for NGG 7321 are shown in Figure |6] where the 
blue curve is Vt, the red curve is 14,t, and the green curve 
is 14,r. The disk geometry parameters agree for both the 
rotation only and bisymmetric models (see Table [HI. This is 
one example where forcing the bisymmetric velocity compo¬ 
nents ( 14,4 and 14,r) to zero for r > 20” improved the fit by 
reducing uncertainties in the resulting velocities. Given the 
lack of coherence in the model residuals, it was determined 
that the bisymmetric model is optimal for this galaxy. This 
is corroborated by X^fC > for the bisymmetric flows in 
Figure [6l 

The bisymmetric model inclination found by DiskFit is 
ife = 46 ± 2° which agrees well with the photometric value 


© 2015 RAS, MNRAS 000, [T]-?? 



















































The incidence of bar-like kinematic flows in CALIFA galaxies 9 



40 20 0 -20 40 40 20 0 -20 -40 

Aa (arcsec) Aa (arcsec) 


Figure 3. DiskFit model results for NGC 1056, a galaxy that is representative of the Rotation Only category, (a) Ho velocity field, (b) 
Uncertainties in the velocity, (c) Best fitting rotation only model, (d) Rotation only model residuals, (e) Best fitting bisymmetric model, 
(f) Bisymmetric model residuals. In all panels, the colourbar is in hms"*^. 


derived within the CALIFA collaboration of ip = 47 ± 1° 
(see Table [P and is higher than the value from 2MASS 
of ip = 39° JSkrutskie et al.ll2006l ~l. The PA, = 12 ± 1° 
is lower than the value of ^ = 25° estimated from 
2MASS data llSkrutskie et al.l l2006ll but agrees well with 


other kinematic works such as 


Barrer a-Ballesteros et al.l 
I 2 OI 4 I) , = and 
= 14.4 ± 2.8° respectively). The morphological PA how- 


ll2014h and iGarcfa-Lorenzo et al 


ever, can be easily shifted by the edges of the spiral 
arms, which lie outside the CALIFA FOV. The systemic 
velocity found is Fsps = 7123 ± 3 kms”*^ (Table [3]) which 
is c omparable, although not c onsistent with, the value 
from iGiovanelli fc Havnesl (Il993l l of Vs vs = 7145 ± 5kms~^ 

i Table HI) . This could be because the (Giovanelli fc Havnej 
I 993 II Vays stems from single-dish HI observations rather 
than detailed modelling of the velocity field: one therefore 
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Figure 5. DiskFit model results for NGC 7321, a galaxy that is representative of the Non-Circular Flows class, (a) Ha velocity field, 
(b) Uncertainties in the velocity, (c) Best fitting rotation only model, (d) Rotation only model residuals, (e) Best fitting bisymmetric 
model, (f) Bisymmetric model residuals. In all panels, the colourbar is in kms“^. 


expects the DiskFit value to be more accurate. The re¬ 
sulting value for the PA of the bar was determined to be 
(/)(, = 48 ± 6° in the sky plane, which seems plausible given 
the orientation of the optical bar in Figure [TJ:. 


4.3 Sample Results 

We use the method outlined in 1)4. H and demonstrated in 1)4.21 
to identify the optimal DiskFit model for each of the 37/100 
CALIFA DRl galaxies that meet the selection criteria of ^ 
for which valid models were obtained. We henceforth restrict 
our analysis to these 37 systems, and discuss the modelling 
results from this final sample as a whole in this section. 

The model classifications in Table [2] show that the H a 
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Figure 6. Kinematic components for NGC 7321, whose optimal 
kinematic model is the bisymmetric model. The blue curve is Vt 
(the rotation curve) for the bisymmetric model shown in Figure 
[S] the red curve is Vi.t a-nd green curve is V 2 ,r- 

velocity fields of 25/37 (67.6j(g'5%) galaxies are best char¬ 
acterized as Rotation Only. Of these, VIl‘2'o (68.0j(^(/4%) 
had an optimal model of rotation only, whereas 8/25 
galaxies had insufficient Ha or suffered from 
spatial masking near the galaxy centre that prevented a ro¬ 
bust search for non-circular flows and were placed in the 
Can’t Tell subcategory. The remaining Vlj'i'l (32.4lgg%) 
galaxies were found to contain Non-Circular Flows, such 
that 10 galaxies were best fit with a bisymmetric model, 
one galaxy with a warped disk model (NGC 36) and one 
galaxy with a radial model (UGC 9476). We address these 
special cases in 115.11 

Figure [7] compares the disk geometry obtained from the 
rotation only and bisymmetric DiskFit models of the same 
galaxy. In all panels blue symbols represent galaxies in the 
Rotation Only category including those classified as Can’t 
Tell, while red symbols show galaxies in the Non-Circular 
Flows category. There are two galaxies for which significant 
discrepancies exist between the best fitting rotation only and 
bisymmetric model disk geometries. First, the linear feature 
1" west of the optical centre in UGC 8733 (see Figure A16) 
was fit differently by various DiskFit models, producing out¬ 
liers in Figure [71 Second, the kinematic inclination of the 
bisymmetric model for UGC 11717 corresponds to the min¬ 
imum value allowed by DiskFit, implying an unreliable fit 
for this galaxy (note that the rotation only model for this 
galaxy is reliable). Figure [7] shows that in general, the disk 
geometry of a given system is consistently determined for all 
the different models: in most cases, it is the amplitude of the 
non-circular flows and residual pattern that determines the 
optimal model, not the disk geometry. This was explicitly 
shown for the case studies in 44.21 This result is also broadly 
consis tent with the conclusions of iBarrera-Ballesteros et al.l 
ll2014l) . who find no significant deviation in the PA of the line 
of nodes in barred galaxies relative to unbarred systems. 

Figure [8] compares the disk geometry of the optimal 
DiskFit model for each galaxy to the photometric litera¬ 
ture values from Table [T] The color coding in Figure [H]is the 
same as that in Figure [T] With the exception of UGC 8733 
discussed above whose difference in the relative kinematic 


Figure 9. Mean value of V 2 ,t (stars) and V 2 , 7 - (triangles) as a 
function of the number of independent radial bins (bottom) or 
arcsec from the centre (top) over which the flow was detected, in 
final sample galaxies where the bisymmetric model was judged 
to be optimal (col (3) = ‘m2’ in Table[2]l. The error bars are the 
standard deviation on the mean. At least one component of each 
galaxy lies above and to the right of the horizontal and vertical 
black dashed lines at U ~ 15kms“'^ and 2 radial bins: we adopt 
these values as the detectability threshold for bisymmetric flows 
in the final sample. 


centre lies beyond the limits of the plot, all of the scat¬ 
ter lies within a single 2.7”-diameter PPaK fibre. There is 
good agreement between the kinematic and photometric PA 
(Figure (Sjt). In a few cases the kinematic PA differs by 20° 
from the photometric PA, however the corresponding disks 
are at the low end of the sample inclination range. Figure [SJ: 
shows a scatter of the photometric inclination derived by the 
CALIFA collaboration compared to the kinematic inclina¬ 
tion from DiskFit. There is good agreement between the val¬ 
ues within uncertainties with the exception of NGG 6154 and 
NGC 7025 (Figures A24 and A32, respectively), for which 
the kinematic inclination is much higher than that returned 
from the photometry. In general, there is good agreement 
between the literature and DiskFit Vaya- There are some 
cases where the DiskFit Vaya differs by ~ 100 kms“^ from 
the literature value, however the latter are highly uncertain 
in these cases (see Table [1]). 

Figure [9] shows the weighted mean value of V 2 ,t (stars) 
and V 2 ,r (triangles) as a function of the number of in¬ 
dependent radial bins (bottom) or arcsec from the cen¬ 
tre (top) over which the non-circular flow was detected for 
those final sample galaxies where the bisymmetric model 
was deemed optimal. The mean value of Vi.t or Vi.r exceeds 
15 kms“^ over two independent radial rings in each galaxy, 
as illustrated by the horizontal and vertical dashed lines 
in Figure [9] We adopt this pair of values as DiskFit’s de¬ 
tectability threshold for non-circular flows in the final sam¬ 
ple, and discuss its implications in 45.21 

Figure [TO] shows the distribution of modelled galaxies as 
a function of morphological type (left) and photometric bar 
classification (right). It illustrates that the photometric bar 
classification agrees with the model selected in 86.2(tg g% 
(25/29) of the galaxies that contained enough Ha to search 
for non-circular flows. Of the 16/37 galaxies photometrically 
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Figure 7. Comparison of DiskFit disk geometry for best-fitting rotation only vs non-circular flow models for each final sample galaxy. 
In all panels, blue dots represent galaxies in the Rotation Only category including those classified as Can’t Tell^ while red dots represent 
galaxies in the Non-Circular Flows category, (a) Position r = \/x'^ +2/^ of the kinematic centre relative to centre pixel, (b) PA. The 
uncertainties are smaller than the symbol size for many systems, (c) Kinematic inclination, ij^. (d) Difference between Vsys returned 
from from the rotation only and bisymmetric models, as a function of that returned by the rotation only model. In all panels, the solid 
black line shows the 1:1 relation between parameters. Outliers from this relation are labelled. 
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Figure 8. Comparison of optimal DiskFit disk geometry with literature values from Tabled for each final sample galaxy. In all panels, 
blue dots represent galaxies in the Rotation Only category including those classified as Can’t Tell, while red dots represent galaxies in 
the Non-Circular Flows category, (a) Difference in centre position of right ascension and declination relative to photometric centre, (b) 
PA. The uncertainties are smaller than the symbol size for many systems, (d) Kinematic inclination (i^) compared with the photometric 
inclination derived by the CALIFA collaboration (ip), (d) Difference between Vsys returned from from the rotation only and bisymmetric 
models, as a function of literature values. In all panels, the solid black line shows the 1:1 relation between parameters. 
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Figure 10. Optimal kinematic model for final sample galaxies as a function of morphological type (left) and bar classification (right). 
The dark blue colour represents the Rotation Only models, light blue are the galaxies classified as Can’t Tell, dark green are cases of 
Non-Circular Flows: bisymmetric, radial and warp models combined. 
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classified as unbarred, 81.2l®3^j% (13/16) were found to be 
best fit with the rotation only model. Of the 14/37 galaxies 
that were photometrically classified as barred, 57 . 11^3 3 % 
(8/14) were found to have significant non-circular flows and 
an additional 35 . 7 l^g' 2 % (5/14) of these galaxies were placed 
into the Can’t Tell subcategory due to the lack of Ha or 
spatial masking near the galaxy centre. Thus while there is 
a good correspondence between the photometric and kine¬ 
matic bar classifications in the final sample, there are excep¬ 
tions. We discuss these statistics in light of the kinematic bar 
detectability implied by Figure [9] in 115.31 


5 DISCUSSION 

In this section we interpret the DiskFit results for the final 
sample of 37/100 DRl galaxies presented in >14.31 A dis¬ 
cussion of the origin of detected non-circular flows is pre¬ 
sented in gsn DiskFit’s sensitivity to non-circular flows is 
addressed in ^5.21 and a comparison of the photometric and 
kinematic classifications of the final sample galaxies is given 
in >15.31 


5.1 Origin of Detected Non-Circular Flows 

By definition (Eq.[2]), DiskFit’s bisymmetric model is sensi¬ 
tive to a bisymmetric flow with constant phase. As discussed 
in >J^ this model is therefore well-suited to detecting bar-like 
flows in disk galaxies. Figure [ 10 ] shows that we detect bisym¬ 
metric flows in most of the photometrically barred galaxies 
in the final sample that have sufficient H a to afford a search. 
Moreover, a visual inspection shows that the kinematic bar 
angle returned by DiskFit in those galaxies (Tabled) as well 
as the extent of the non-circular flow regions agree well with 
the photometric bar properties observed in the SDSS igr 
images. We therefore conclude that the bisymmetric flows 
detected by DiskFit do indeed represent physical bar-like 
flows, and interpret the main results in >|4]in this context. 

There are, however, 2/37 (5.4jl® g%) sample galaxies for 
which another non-circular flow model is preferred to the 
bisymmetric model: a radial flow model was deemed optimal 
in UGC 9476 (Fig. A20), whereas a warped disk was invoked 
for NGC 36 (Fig. A3). We discuss these special cases here. 

Pure radial flows ha ve not been observed (e.g 
IWong. Blitz fc Bosmal l2004ll and are not physically well- 
motivated. It seems unlikely that the radial flows detected 
in UGG 9476 are indeed pure radial flows, but that Disk- 
Fit has rather mis-identified a bar-like flow. This is possi¬ 
ble when the bar is aligned along either the major or mi¬ 
nor axis of the disk {9h —>■ 0° and 6b —^ 90° in Eq. [2|), 
since the non-circular flow components in the bisymmet¬ 
ric model b ecome degenerate with rot ation at these bar 
orientations (ISellwood &: Sanched 1201011 . The radial model 
(Eq. (31 does not have this same degeneracy and is more re¬ 
liable in these circumstances (recall that both m = 0 and 
m = 2 flows in the disk plane project to m = 1 flow s 
in the sky plane; ISchoenmakers. Franx fc de Zeeuwill997ll . 
Utilizing the radial model in cases where the bar angle ap¬ 
proaches the major or minor axis allows for the continued 
detection of non-circular flows for these bar geometries. We 
propose that there is a kinematic bar signature in UGG 9476 
along its major axis, which has been detected by the radial 


model. Despite its classification as photometrically unbarred 
(Table [T]), a linear feature appears to be visible along the 
major axis in the SDSS igr composite image (see Fig. A20) 
that approximately matches the extent of the radial flows 
detected. A more detailed analysis of this galaxy will be per¬ 
formed using DiskFit’s photometric branch as part of future 
work. 

We find that the kinematics of NGC 36 are best de¬ 
scribed by a warped disk model instead of non-circular flows 
in a flat disk. This galaxy is photometrically classified as 
barred however, we find significant uncertainties in the de¬ 
rived velocities for the bisymmetric model. There are hints 
of a radial change in PA in the velocity field thus indicating 
that a warp may be present. It is also possible that NGC 36 
hosts an oval disk instead of a warp, which woul d have a sim¬ 
ilar kinematic signature (e.g. lKormendvll2013l '). If NGC 36 
has an inner warp, the frequency of 1/37 (2.7 ln'a%) galax¬ 
ies seems high compared to other studies fe.g. [Brigg sl 199d; 
iGarcfa-Ruiz. Sancisi fc KuiikerJ I2OO2I : Ide Blok et al.il2008l j: 
however, none of these studies selected statistically repre¬ 
sentative samples. 

The discussion above suggests that the warped flow that 
we detect in NGC 36 is not a mis-identified bar-like flow, but 
results from different physical processes in the disk. This 
suggests that 1/12 ( 8 . 312 ^ 3 ®%) ^^e coherent non-circular 
flows in intermediate-late type CALIFA galaxies do not stem 
from bars, with the important caveat that the sample stud¬ 
ied here is relatively small. Since our discussion in the sec¬ 
tions below pertains to bar-like non-circular flows, we omit 
NGC 36 from further consideration. 


5.2 Sensitivity to Bar-Like Flows 

Figure ini illustrates that the detected non-circular flows in 
the final sample have at least one component with an ampli¬ 
tude that exceeds 15kms“^ over at least two independent 
radial bins. As discussed in >15.11 the physical implication of 
this threshold for bar-like flows in galaxies depends on the 
bar angle: Eq. [3] becomes degenerate for bars close to either 
the major or minor axis of the disk, making them more dif¬ 
ficult to reliably detect. A thorough simulation of the effect 
for the CALIFA sample is beyond the scope of this work, 
and isn’t justified given the relatively small sample studied 
here. Our experience suggests, however, that DiskFit’s sen¬ 
sitivity to bar-like flows is unaffected by bar angle when the 
latter is more than ~ 10 ° from the disk major or minor axis 
in intermediate-inclination galaxies. The bar-like flows de¬ 
tected by DiskFit for the 10 galaxies in Figure [Sjare at least 
20 ° from either the major or minor axis, suggesting that 
bar angle is not influencing the flow detectability therein. 

Of the 8/37 ( 21 . 6 jlg' 2 %) gas-rich galaxies photomet¬ 
rically classified as barred for which the bisymmetric 
model is optimal, the average detected flow amplitude is 
V2,t,avg = 9kms“^ Over ~ 3.5kpc and V2,r,avg = 28kms“^ 
over 5.5 kpc. Small-number statistics preclude us from 
corre cting those number s for the CALIFA selection func¬ 
tion (IWalcher et al.l^2014^ . However, the fact that we detect 
such flows in most (8/9 or 88 . 91 ^ 3 ^ 3 %) of the photometri¬ 
cally barred galaxies with sufficient Ha to afford a search 
suggests that this flow amplitude is characteristic of barred 
galaxies with the masses probed by CALIFA. We conclude 
that non-circular flows due to bars in CALIFA galaxies are 
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detectable by DiskFit, provided that they are at intermedi¬ 
ate angles to the disk major and minor axes. 

On the other hand, we detect non-circular flows in only 
1/7 ( 14 . 345 ^ 3 ^%) galaxies of photometric classihcation AB. 
This galaxy (NGC 477; Figure A4) has non-circular flows 
as high as 70kms“^, extending ~ 4kpc from the center. 
There is a masked star near the minor axis of the velocity 
field that could have interfered with the fit, raising the pos¬ 
sibility that the non-circular flows in this system are over¬ 
estimated. Regardless, our non-detection of most AB galax¬ 
ies suggests that the non-circular flows associated with their 
weaker/intermediate bars fall below our detection threshold. 
At the characteristic distance of these galaxies in our sam¬ 
ple, this implies that these flows are weaker than 15kms“^ 
and/or do not extend more than 2.25 kpc in radius. 

5.3 Galaxies with Different Kinematic and 
Photometric Bar Classifications 

Figure [To] illustrates that for the majority of the galaxies in 
the final sample, the kinematic classification obtained from 
our DiskFit models matches the photometric classification 
in Table [T] By and large, we find bar-like flows in photomet¬ 
rically barred galaxies with sufficient H a to afford a search, 
and that the rotation only model best describes the veloc¬ 
ity fields of photometrically unbarred galaxies. There are 
exceptions, however: we find no evidence of bar-like flows 
in the barred galaxy UGC 5359 (Fig. A9), and signihcant 
non-circular flows in the photometrically unbarred galaxies 
NGG 7819, UGG 9476, and NGG 7025 (Figs. Al, A20, and 
A32, respectively). In this section, we investigate the origin 
of these discrepancies to estimate the incidence of systems 
where the photometric bar classification is different from 
that inferred kinematically. 

We find no evidence for bar-like flows in the photomet¬ 
rically barred galaxy UGG 5359 (Fig. A9). However, an ex¬ 
amination of the SDSS image reveals that the bar angle 
approaches the disk major axis, making bar-like flows dif¬ 
ficult to detect (see >15.11 and >15.211 . It is therefore plausible 
that the amplitude of the bar-like flows in UGC 5359 resem¬ 
ble those of other barred galaxies in the sample, but that 
the bar angle precludes detecting these flows with DiskFit’s 
bisymmetric model. We see some evidence for non-circular 
flows in the radial model for this galaxy, although they are 
not as convincing as those for UGC 9476 in §§unusual2. 

On the other hand, we do detect non-circular 
flows in the photometrically unbarred galaxies NGC 7819, 
UGC 9476, and NGC 7025 (Figs. Al, A20, and A32, respec¬ 
tively). We discussed UGC 9476 in >)5.1l suggesting that it 
contains a major axis bar and is therefore photometrically 
mis-classified. For NGC 7819 and NGC 7025, the average 
non-circular flows detected are V 2 ,t,avg = 85 kms“^ over ~ 8 
3"rings (~ 9 kpc) and V 2 ,r,avg = 48kms“^ over ~ 6 3"rings 
(~ 6.5 kpc): Figure (HI shows that they are well above our 
detection threshold. We re-examined the SDSS images for 
these two systems and find no evidence for a bar that was 
missed during the photometric classihcation. It is possible 
that the non-circular hows that we detect are being driven 
by anot her mechani s m, su ch as an interaction with a nearby 
galaxy. iMahtessianI lll998ll classify NGC 7819 as a member 
of a group by searching for nei ghbours with similar radial 
velocities. NGC 7025 is isolated jKarachentsevairi973l l . how¬ 


ever, and we see no obvious photometric feature that could 
drive the non-circular hows that we detect. It is possible 
that a bar-like feature in NGC 7025 is present in the NIR but 
not in the optical. This possibility is unlikely, however, since 
the amplitude of the non-circular hows implies a relatively 
strong bar, and weak b ars are the ones that are typically 
obscured in the optical llMarinova et al.ll2009al l. 

Considering the above discussion, and excluding inter¬ 
mediate (AB) bars, we conservatively conclude that the pho¬ 
tometric classihcation of the hnal sample galaxies with suf- 
hcient H a to enable a search for non-circular hows belies a 
different kinematic classihcation in at least 2/23 (8.7l2’g%) 
systems (UGC 9476 and NGC 7025). Systematic searches 
for bar-like hows thus not only characterise the kinematic 
properties of galaxy bars, but may also reveal galaxies in 
which bar-like hows are driving galaxy evolution despite the 
lack of a clear photometric bar signature. 

Our analysis of the 37/100 CALIFA DRl galaxies suit¬ 
able for kinematic modelling hints at the richness of the 
information that can be gleaned from a systematic search 
for non-circular hows in nearby galaxies, but the relatively 
small hnal sample size precludes a detailed statistical inter¬ 
pretation of our results. That is set to change with the full 
CALIFA sample of ~ 600 galaxies that will be available in 
the near future: scaling the results presented here, we expect 
that ~ 200 CALIFA galaxies will ultimately afford detailed 
kinematic analyses. A joint photometric and kinematic de¬ 
composition of the SDSS images for each of these galaxies 
would enable quantitative comparisons between the bars de¬ 
tected in galaxy images and the non-circular hows found in 
their velocity helds, affording a three-dimensional examina¬ 
tion of nearby barred systems. This work is underway. 


6 SUMMARY 

We have used DiskFit to model the H a velocity helds of 
gas-rich, intermediate-inclination CALIFA DRl disk galax¬ 
ies in a direct search for bar-like hows. We apply rotation 
only, bisymmetric hows, and (in some cases) radial how and 
warped disk models to 49/100 galaxies with photometric in¬ 
clinations 40° < ip < 70° and a visible velocity held gradi¬ 
ent, and hnd acceptable models for a hnal sample of 37/100 
systems. For each galaxy in the hnal sample, we use a 
test to search for statistically signihcant non-circular hows 
as well as examine residual plots to determine the optimal 
kinematic model. 

We hnd good agreement between the disk geometry re¬ 
turned by DiskFit for diherent models of the same galaxy, 
as well as between the optimal model values and the pho¬ 
tometric disk geometries and systemic velocities from the 
literature. Of the 29/37 hnal sample systems with sufficient 
H a near the galaxy centre to afford a search for non-circular 
hows, we deem that 17/29 (58.61®' 4 %) are best described by 
the rotation only model while 12/29 (41.4lg'3%) contain sta¬ 
tistically signihcant non-circular hows. Of these latter galax¬ 
ies we hnd that the bisymmetric model is optimal for 10/12 
(83.3l45®5%) systems, and favour radial hows and a warped 
disk for the remaining two galaxies, respectively. At least 
one bisymmetric how component in each of the 10 galax¬ 
ies exceeds 15kms“^ over at least two independent radial 
bins (~ 2.25 kpc at the characteristic hnal sample distance of 
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77Mpc): we adopt this pair of values as the detectability 
threshold for DiskFit in the final sample. 

Accounting for the low sensitivity of the bisymmetric 
model to bar-like flows aligned near the major or minor axis 
of the disk and comparing our optimal kinematic models to 
the photometric bar classifications of the sample galaxies, 
we conclude that the non-circular flows that we detect in 
11/12 (91.7lj4®g%) galaxies stem from bars, while the re¬ 
maining system likely harbours an inner warp or oval disk. 
We find that photometrically barred CALIFA DRl galaxies 
have an average non-circular flow (Vi.t and Vi.r) amplitude 
of Vi = 16 ± lkms“^ over a radial extent of 4.5 kpc. On 
the other hand, the absence of non-circular flows in galaxies 
with intermediate (AB) bar classifications implies that these 
flows fall below our detection threshold for intermediate bar 
angles. It is evident from these relations that bisymmetric 
flows are nearly ubiquitous in strongly barred galaxies, as 
one would expect. 

There are 4 galaxies in the final sample where our kine¬ 
matic classification differs from that obtained photometri¬ 
cally: we find that the barred galaxy UGC 5359 is best de¬ 
scribed by a rotation only model, and detect bar-like non¬ 
circular flows in the unbarred systems NGC 7819, UGC 9476 
and NGC 7025. We find it plausible that DiskFit missed the 
putative bar-like flows in UGC 5359 due to geometric effects, 
and that the photometric classification of UGC 9476 is in¬ 
correct. It is also plausible that interactions in the group in 
which NGC 7819 resides have caused the non-circular flows 
that we detect; however, NGC 7025 is isolated. We there¬ 
fore conclude that in 2/23 (8.7/l2'g%) galaxies - or ~ 10% 
of the time - the photometric classification of CALIFA DRl 
galaxies belies a different kinematic classification. 


ACKNOWLEDGEMENTS 

KS acknowledges support from the Natural Sciences and En¬ 
gineering Research Council of Canada. CJW acknowledges 
support through the Marie Curie Career Integration Grant 
303912. RAM is funded by the Spanish program of Interna¬ 
tional Campus of Excellence Moncloa (CEI). JMA acknowl¬ 
edges support from the European Research Council Starting 
Grant (SEDmorph; P.I. V. Wild). 

This study makes uses of the data provided by the 
Calar Alto Legacy Integral Field Area (CALIFA) sur¬ 
vey (http://califa.caha.es/ ). Based on observations col¬ 
lected at the Centro Astronomico Hispano Aleman (CAHA) 
at Calar Alto, operated jointly by the Max-Planck-Institut 
fur Astronomie and the Institute de Astrofisica de Andalu- 
cia (CSIC). 

This research has made use of the NASA/IPAC Ex- 
tragalactic Database (NED) which is operated by the Jet 
Propulsion Laboratory, California Institute of Technology, 
under contract with the National Aeronautics and Space 
Administration. 

Funding for the SDSS and SDSS-II has been provided 
by the Alfred P. Sloan Foundation, the Participating In¬ 
stitutions, the National Science Foundation, the U.S. De¬ 
partment of Energy, the National Aeronautics and Space 
Administration, the Japanese Monbukagakusho, the Max 
Planck Society, and the Higher Education Funding Council 
for England. The SDSS Web Site is http: //www. sdss. org/ 


The SDSS is managed by the Astrophysical Research 
Consortium for the Participating Institutions. The Partic¬ 
ipating Institutions are the American Museum of Natu¬ 
ral History, Astrophysical Institute Potsdam, University of 
Basel, University of Cambridge, Case Western Reserve Uni¬ 
versity, University of Chicago, Drexel University, Fermilab, 
the Institute for Advanced Study, the Japan Participation 
Group, Johns Hopkins University, the Joint Institute for 
Nuclear Astrophysics, the Kavli Institute for Particle As¬ 
trophysics and Cosmology, the Korean Scientist Group, the 
Chinese Academy of Sciences (LAMOST), Los Alamos Na¬ 
tional Laboratory, the Max-Planck-Institute for Astronomy 
(MPIA), the Max-Planck-Institute for Astrophysics (MPA), 
New Mexico State University, Ohio State University, Uni¬ 
versity of Pittsburgh, University of Portsmouth, Princeton 
University, the United States Naval Observatory, and the 
University of Washington. 


REFERENCES 

Abazajian K. et ah, 2009, ApJS, 182, 543 
Abazajian K. et ah, 2003, AJ, 126, 2081 
Adelman-McCarthy J. K. et ah, 2008, ApJS, 175, 297 
Aguerri J. A. L., Mendez-Abreu J., Corsini E. M., 2009, 
A&A, 495, 491 

Aguerri J. A. L. et ah, 2015, preprint (arXiv:1501.05498) 
Athanassoula E., 1992, MNRAS, 259, 345 
Athanassoula E., Machado R. E. G., Rodionov S. A., 2013, 
MNRAS, 429, 1949 

Barazza F. D., Jogee S., Marinova L, 2008, ApJ, 675, 1194 
Barrera-Ballesteros J. K. et ah, 2014, A&A, 568, A70 
Bershady M. A., Verheijen M. A. W., Swaters R. A., An¬ 
dersen D. R., Westfall K. B., Martinsson T., 2010, ApJ, 
716, 198 

Blanc G. A. et al., 2013, AJ, 145, 138 
Briggs F. H., 1990, ApJ, 352, 15 
Bryant J. J. et al., 2015, MNRAS, 447, 2857 
Bundy K. et al., 2015, ApJ, 798, 7 
Cameron E., 2011, PASA, 28, 128 
Cappellari M. et al., 2011, MNRAS, 413, 813 
Croom S. M. et al., 2012, MNRAS, 421, 872 
Davoust E., Considere S., 1995, A&AS, 110, 19 
de Blok W. J. G., Walter F., Brinks E., Trachternach C., 
Oh S.-H., Kennicutt, Jr. R. C., 2008, AJ, 136, 2648 
de Vaucouleurs G., de Vaucouleurs A., Corwin J. H. G., 
Buta R. J., Paturel G., Fouque P., 1991, Third Reference 
Catalogue of Bright Galaxies, Version 3.9 
de Zeeuw P. T. et al., 2002, MNRAS, 329, 513 
deirAntonio 1. P., Bothun G. D., Geller M. J., 1996, AJ, 
112, 1759 

Epinat B., Amram P., Marcelin M., 2008, MNRAS, 390, 
466 

Epinat B. et al., 2008, MNRAS, 388, 500 
Erwin P., Debattista V. P., 2013, MNRAS, 431, 3060 
Eskridge P. B. et al., 2000, AJ, 119, 536 
Falco E. E. et al., 1999, PASP, 111, 438 
Freudling W., Haynes M. P., Giovanelli R., 1992, ApJS, 79, 
157 

Friedli D., Martinet L., 1993, A&A, 277, 27 
Garci'a-Lorenzo B. et al., 2014, A&A, 573, A59 


© 2015 RAS, MNRAS 000, [T]-?? 


16 L. Holmes et al. 


Garcfa-Ruiz I., Sancisi R., Kuijken K., 2002, A&A, 394, 
769 

Giovanelli R., Haynes M. P., 1993, AJ, 105, 1271 
Haynes M. P., Giovanelli R., Herter T., Vogt N. P., 
Freudling W., Maia M. A. G., Salzer J. J., Wegner G., 
1997, AJ, 113, 1197 

Hernandez O., Carignan G., Amram P., Ghemin L., Daigle 
O., 2005, MNRAS, 360, 1201 

Hnchra J. P., Geller M. J., Corwin, Jr. H. G., 1995, ApJS, 
99, 391 

Hnsemann B. et al., 2013, A&A, 549, A87 
Karachentseva V. E., 1973, Astrofizicheskie Issledovaniia 
Izvestiya Spetsial’noj Astrofizicheskoj Observatorii, 8, 3 
Kehrig C. et al., 2012, A&A, 540, All 
Kelz A. et al., 2006, PASP, 118, 129 
Kent B. R. et al., 2008, AJ, 136, 713 
Kormendy J., 2013, Secular Evolution in Disk Galaxies, 
Falcon-Barroso J., Knapen J. H., eds., p. 1 
Kormendy J., Bender R., 2012, ApJS, 198, 2 
Kormendy J., Kennicutt R. C., 2004, ARA&A, 42, 603683 
Kuzio de Naray R., Arsenault C. A., Spekkens K., Sell- 
wood J. A., McDonald M., Simon J. D., Teuben P., 2012, 
MNRAS, 427, 2523 

Laurikainen E., Salo H., Athanassoula E., Bosma A., 
Herrera-Endoqui M., 2014, MNRAS, 444, L80 
Lu N. Y., Hoffman G. L., Groff T., Roos T., Lamphier C., 
1993, ApJS, 88, 383 
Mahtessian A. P., 1998, Afz, 41, 308 

Marinova I. et al., 2009a, in Astronomical Society of the 
Pacific Conference Series, Vol. 419, Galaxy Evolution: 
Emerging Insights and Future Challenges, Jogee S., Mari¬ 
nova I., Hao L., Blanc G. A., eds., p. 138 
Marinova I. et al., 2009b, ApJ, 698, 1639 
Marzke R. O., Hnchra J. P., Geller M. J., 1996, AJ, 112, 
1803 

Masters K. L. et al., 2011, MNRAS, 411, 2026 
Mendez-Abreu J., Debattista V. P., Corsini E. M., Aguerri 
J. A. L., 2014, A&A, 572, A25 
Mulchaey J. S., Regan M. W., 1997, ApJL, 482, L135 
Nishiura S., Shimada M., Ohyama Y., Murayama T., 
Taniguchi Y., 2000, AJ, 120, 1691 
Nordgren T. E., Chengalur J. N., Salpeter E. E., Terzian 
Y., 1997, AJ, 114, 77 
O’Neil K., 2004, AJ, 128, 2080 

Ramella M., Geller M. J., Hnchra J. P., Thorstensen J. R., 
1995, AJ, 109, 1458 

Reese A. S., Williams T. B., Sellwood J. A., Barnes E. I., 
Powell B. A., 2007, AJ, 133, 2846 
Regan M. W., Sheth K., Vogel S. N., 1999, ApJ, 526, 97 
Rosales-Ortega F. F., Kennicutt R. C., Sanchez S. F., Dfaz 
A. I., Pasquali A., Johnson B. D., Hao C. N., 2010, MN¬ 
RAS, 405, 735 

Roth M. M. et al., 2005, PASP, 117, 620 
Saintonge A., Giovanelli R., Haynes M. P., Hoffman G. L., 
Kent B. R., Martin A. M., Stierwalt S., Brosch N., 2008, 
AJ, 135, 588 

Sanchez S. F., 2006, Astron. Nachrichten, 327, 850 
Sanchez S. F. et al., 2012, A&A, 538, A8 
Schoenmakers R. H. M., Franx M., de Zeeuw P. T., 1997, 
MNRAS, 292, 349 

Sellwood J. A., Sanchez R. Z., 2010, MNRAS, 404, 1733 
Sellwood J. A., Wilkinson A., 1993, Reports on Progress 


in Physics, 56, 173 

Sheth K., Mizusawa T., Kim T., Munoz-Mateos J., Regan 
M. W., de Swardt B., Gadotti D., S4G Team, 2011, in 
Bulletin of the AAS, Vol. 43, A AS Meeting Abstracts # 
217, p. 207.06 

Sheth K., Vogel S. N., Regan M. W., Thornley M. D., 
Teuben P. J., 2005, ApJ, 632, 217 
Shlosman L, Begelman M. C., Frank J., 1990, Nature, 345, 
679 

Shlosman L, Frank J., Begelman M. C., 1989, Nature, 338, 
45 

Simard L., Mendel J. T., Patton D. R., Ellison S. L., Mc- 
Connachie A. W., 2011, ApJS, 196, 11 
Skrutskie M. F. et al., 2006, AJ, 131, 1163 
Spekkens K., Sellwood J. A., 2007, ApJ, 664, 204 
Theureau G., Bottinelli L., Coudreau-Durand N., 
Gouguenheim L., Hallet N., Loulergue M., Paturel G., 
Teerikorpi P., 1998, A&AS, 130, 333 
Theureau G., Hanski M. O., Coudreau N., Hallet N., Mar¬ 
tin J.-M., 2007, A&A, 465, 71 
Walcher C. J. et al., 2014, A&A, 569, Al 
Wegner G. et al., 2003, AJ, 126, 2268 
Wegner G., Haynes M. P., Giovanelli R., 1993, AJ, 105, 
1251 

Weinberg M. D., Katz N., 2002, ApJ, 580, 627 
Weiner B. J., Sellwood J. A., Williams T. B., 2001, ApJ, 
546, 931 

Whyte L. F., Abraham R. G., Merrifield M. R., Eskridge 
P. B., Frogel J. A., Pogge R. W., 2002, MNRAS, 336, 1281 
Wong T., Blitz L., Bosma A., 2004, ApJ, 605, 183 


© 2015 RAS, MNRAS 000, [T]-?? 


